Fidelity, dynamic structure factor, and susceptibility in critical phenomena 
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Motivated by the growing importance of fidelity in quantum critical phenomena, we establish a 
general relation between fidelity and structure factor of the driving term in a Hamiltonian through 
a newly introduced concept: fidelity susceptibility. Our discovery, as shown by some examples, 
facilitates the evaluation of fidelity in terms of susceptibility using well developed techniques such 
as density matrix renormalization group for the ground state, or Monte Carlo simulations for the 
states in thermal equilibrium. 
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I. INTRODUCTION 

Recently, much attention P, [H, H, IH H, @ has been 
drawn to the role of fidelity a concept emerged from 
quantum information theory [7[ , in quantum critical phe- 
nomena Q . Since fidelity is a measure of similarity be- 
tween states, a dramatic change in the structure of the 
ground state around quantum critical point should re- 
sult in a great difference between the two ground states 
on the both sides of critical point. For example, in one- 
dimensional XY model, fidelity shows a narrow trough at 
phase transition point [2|. Similar properties were also 
found in fermionic [|| and bosonic systems pjj. As fi- 
delity is purely a quantum information concept, these 
works actually established a connection between quan- 
tum information theory and condensed matter physics. 

However, except a few specific models such as the one- 
dimensional XY model and Dicke modelQ,Q, it is tedious 
to evaluate fidelity from the ground state wavefunctions. 
Therefore, a neater and simpler formalism is of great im- 
portance for the extensive application of fidelity to criti- 
cal phenomena. For this purpose, we introduce the con- 
cept of fidelity susceptibility, which defines the response 
of fidelity to the driving parameter of the Hamiltonian. 
At zero temperature, we show that fidelity susceptibility 
is intrinsically related to the dynamic structure factor of 
the driving Hamiltonian, namely Hj, that causes quan- 
tum phase transition. Based on some well-developed nu- 
merical techniques for the ground state properties such as 
exact diagonalization (ED) M and density matrix renor- 
malization group (DMRG) [Hj, an applicable scheme is 
proposed to evaluate the dynamic structure factor of Hj. 
On the other hand, starting from the definition of fidelity 
of a thermal state, we show that fidelity susceptibility is 
simply the thermal fluctuation term, such as specific heat 
C v for the internal energy and magnetic susceptibility \ 
for magnetization. These can easily be obtained from 
Monte-Carlo simulations 



II. GROUND-STATE FIDELITY AND 
DYNAMIC STRUCTURE FACTOR 

A general Hamiltonian of quantum many-body sys- 
tems reads 



H(X)=H Q + XH I , 



(1) 



where Hj is the driving Hamiltonian and A denotes 
its strength. The eigenstates \^ n (X)) which satisfy 
#(A)|*„(A)} = E„|* n (A)) define a set of orthogonal 
complete basis in the Hilbert space. Here we restrict 
ourselves to the phase transition which is not induced by 
the ground state level-crossing. That means the ground 
state of the Hamiltonian is non-degenerate for a finite 
system. Next we change A — > A + <5A where 5X is so small 
that perturbation is applicable. To the first order, the 
ground state becomes 



|*o(A + 5A)> = |* (A))+<5A^ 



n#0 



H n0 (\)\* n (\)) 
E (X)-E n (X) 



and 



H n0 = (*„(A)|fl>|* (A)). 



(2) 



(3) 



Following Ref. 0, fidelity is defined as the overlap be- 
tween |* (A)) and * (A + c5A)), that is 



F i (A,5) = |<*o(A)|* (A + £))|. 
Therefore, to the lowest order, we have 



Ft 



\-8X 2 



K<MA)|g,|*E (A))|' 
[E n (X) - E (XW 



(4) 



(5) 
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Clearly, fidelity is <5A-dependent and so it is an artificial 
quantity. Despite of this, we can still see from Eq.© that 
the most relevant term in determining fidelity is its sec- 
ond derivative. Compared with linear response theory, 
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the coefficient term before 5X 2 actually defines the re- 
sponse of fidelity to a small change in A. From this point 
of view, we introduce a new concept fidelity susceptibility 
as 



XF 



lim 



-21nK- 



With Eq. ([5]), it can be rewritten as 



(G) 



on the spectra of the Hamiltonian H(X) and the hopping 
matrix H n0 . Unfortunately, except for some very small 
systems which are usually far away from the scaling re- 
gion, Eq.([7]) is almost not computable due to the lack 
of knowledge in the set of eigenstates. In order to over- 
come the difficulty, it is necessary to consider the time 
evolution of the system. For simplicity, we omit the pa- 
rameter A in the following expressions. Define dynamic 
fidelity susceptibility as 



Xf 



|(^ t (A)|g 7 |vl/o(A))| 2 
[E n (X) - So (A)] 2 



(7) 



in the ground state. We would like to point out that 
although the aforementioned procedure is based on per- 
turbation theory, fidelity susceptibility (J7JI only depends 



Xf 



|QE n |ff/|tt )| s 
[E n - E } 2 + u 



(8) 



Perform a Fourier transformation and take derivative, 
we then obtain 



9xf(t) 
dr 



-it [<¥ |iT/(T)fri(Q)|*o) - <*o|ff/|*o) 2 ] 0(t) + n [(* |^/(0)fl>(r)|* > - (*o|^/|*o) 2 ] B(-t), (9) 

I 



with t being the imaginary time and 

HAt) = e H ^ T H ie - H ^ 



(10) 



The two equations mentioned above are impressive as 
they reveal the mystery of fidelity in understanding quan- 
tum critical phenomena. The terms in the brackets in Eq. 
([9]) are nothing but the dynamic structure factor of Hj. 
Therefore in the original definition of fidelity, we subcon- 
sciously choose the driving term Hj as a candidate of the 
order parameter, though we might not think so at that 
time. 

In order to arrive at a more computable formula, we 
carry out an inverse Fourier transformation and obtain 



Xf = J t [<*o|ff/(r)i//(0)|*o) - (*o|.H/|¥o> a ] dr 

(11) 

where the first term in the bracket can be calculated by 



(¥ |#/(t)#/(0)|¥o) 

T n(_-i\n 

= E LA zr LeTE °(*o\H I H n H I \* ). (12) 



n 



Although fidelity is difficult to calculate from the ground 
state wavefunctions, Eq. (fTTj) and (fT2"| provide us with 
another practical way. Especially, Eq. (p~2|) can be easily 
evaluated via the prevailing numerical techniques, say, 
ED and DMRG. For the ED, once the ground state is 
obtained, the map from one state to another new state 



is just a standard Lanczos step. While for the DMRG, 
the standard algorithm involves a transformation of the 
Hamiltonian of a system and its environment, from a set 
of old basis to a set of new basis which is constructed by 
the m largest weighted eigenstates of the reduced density 
matrix. Precisely, for the system block Hl — O^HlOl, 
and environment block H R = O r H r Or, where Ol(H) 
are constructed from m largest weighted eigenstates of 
the corresponding reduced density matrix. The only 
modification is that, in addition to H(X), Hi should be 
independently transformed in the DMRG procedure, i.e. 
~H i.l = 0{H I>L L and Hj, R = ] R H I>R R . Once the 
final ground state is obtained, the mapping |^') = Hj\^) 
and |\&') = H\^f) is simply a standard step. 



III. FIDELITY AND FIDELITY 
SUSCEPTIBILITY IN SOME MODELS 

In previous works [J 0, H, 0, H, [1], the fidelity gives 
us a deep impression that it can be used to signal an 
phase transition. For example, in the Ising model, [l], 
the fidelity collapses to zero at the critical point. From 
Eq. ((HJ), such a singular behavior can be interpreted as 
the divergence of the fidelity susceptibility. This interpre- 
tation is consistent with that the fidelity is intrinsically 
based on the Landau's symmetry-breaking theory. The 
order parameter is just the driving term in the Hamil- 
tonian, then the divergence of the fidelity susceptibility 
becomes a crucial condition to signal a phase transition 
of Landau's type. However, this restriction also make the 
fidelity fails to identify those phase transitions that is of 
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FIG. 1: (Color online) LEFT: the fidelity between two states 
separated by different 8U = 0.04, 0.05, 0.06, versus U/t of the 
half-filled Hubbard model with N = L — 10 and periodic 
boundary conditions. RIGHT: the fidelity susceptibility \f 
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FIG. 2: (Color online) Fidelity susceptibility as a function of 
U/t for various system size. In order to avoid degeneracy in 
the ground state, the systems of L — 6 and L = 10 are diago- 
nalized with periodic boundary conditions, while the systems 
of L = 8 and L — 12 with antiperiodic boundary conditions. 



infinite order. 

To check this point, also verify the above expression, 
we now study the fidelity susceptibility of a non-trivial 
model in condensed matter physics, more specifically, the 
one-dimensional Hubbard model, where the phase tran- 
sition happened at the half-filling case is of infinite order. 
The Hamiltonian of the Hubbard model reads 



H 



E 

Wo 



(13) 



where cj a and Cj )Cr are the creation and annihilation op- 
erators for electrons with spin a (with a =T)|) at site 



j respectively, 



'-j^J.oi 



t is the hopping integral, 
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and U denotes the strength of on-site interaction. At 



FIG. 3: The scaling behavior of the fidelity susceptibility up 
to 1906 sites for the half-filled Hubbard at U = 0. Clearly, 
X/N is a constant if N is very large. 



half filling, the ground state of the Hubbard model un- 
dergoes a quantum phase transition from an ideal con- 
ductor to a Mott-insulator at the point U = [HI]. For 
simplicity, we diagonalize the Hamiltonian for a 10-site 
systems with periodic boundary conditions using Lanczos 
method, and compute fidelities with various interaction 
intervals and their corresponding susceptibility. Numer- 
ical results are shown in Fig. [T] Obviously, they support 
our conclusion that fidelity susceptibility rather than fi- 
delity is more crucial in the ground state. That is the 
fidelity susceptibility does not depend on the value of 
SU, but the fidelity does. This fact makes it possible 
to evaluate the fidelity from the ground state |\E'o(^ r )) 
without knowledge of \^ (U ± 5U/2)). 

Another interesting observation is that neither fidelity 
susceptibility is a maximum, nor fidelity is a minimum at 
the critical point. This fact goes somewhat beyond the 
physical intuition aroused by the original research mo- 
tivation which expects that the fidelity should be a 
minimum at the critical point. Though finite-size scaling 
analysis (see Fig. [5J shows that the fidelity susceptibility 
at U = may become larger and larger as the system size 
increases, the divergence of the rescaled fidelity suscep- 
tibility xf/N [l4[ at the critical point is an unexpected 
phenomenon (See Fig. [3]). ft is because for the half- filled 
Hubbard model, the phase transition at U — is of the 
Kosterlitz-Thouless type. In this case, the ground-state 
energy as a function of U can be expanded to any or- 
der around U = 0, and the density-density correlation 
defined by the U term in the Hamiltonian (fT3"|) does not 
have a long-range order, then the local order parame- 
ter for the symmetry-breaking theory is not well-defined. 
On the other hand, if U deviates from the critical point, 
the scaling behavior shown in Fig. [5] manifests that the 
fidelity susceptibility may increases in the small U re- 
gion. However, as the L increases, the trend of xf is still 
similar to the case of U = 0. Therefore, to find an ap- 
propriate driving term is very crucial for understanding 
the role of fidelity in quantum critical phenomena. For 
the Hubbard model, besides the U term, it seems be very 
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necessary to introduce another non-local term for the fi- 
delity. Nevertheless, it is still a challenging problem at 
present. 



IV. MIXED-STATE FIDELITY AND THERMAL 
PHASE TRANSITIONS 

The generalization of fidelity to finite temperatures is 
proposed recently. Based on the definition of fidelity be- 
tween two mixed states, it has been shown that fidelity 
can be expressed in terms of the partition function Q , 



Z{f3) 



^/Z(P-Sp/2)Z(P + Sp/2Y 
where (3 = 1/T, and 

Z{P) = Y,e-^=Y J 9{E)e-^. 



(14) 



(15) 



Here g(E) is the density of states and can be calculated 
from Monte-Carlo simulations [ll[ such as Wang-Landau 
algorithm [jjj • Then the fidelity susceptibility driven by 
temperature can be calculated as 



-2lnF, 



Xf = 



S[3 2 



8(3^0 



A/3 2 ' 



(16) 



Similarly, if the driving term in the Hamiltonian is 
a Zeemann-like term, which is crucial in Landau's 



symmetry-breaking theory, then the fidelity susceptibil- 
ity is simply the magnetic susceptibility \, 



-21nF, 



XF 



5h 2 



fix 



(17) 



Sh^Q 



Clearly, the specific heat is simply the fluctuation of the 
internal energy, i.e. C v — f3 2 {(E 2 ) - (E) 2 ), while the 
magnetic susceptibility is the fluctuation of the magne- 
tization, i.e. x = /3((M 2 ) - (M) 2 ). Thus fidelity sus- 
ceptibility is just the fluctuation (structure factor) of the 
driving term in the Hamiltonian. 



V. SUMMARY 

In summary, we established a general relation between 
the fidelity and the structure factor of the driving term 
of the Hamiltonian for both quantum and classical crit- 
ical phenomena. Such relation not only enables us to 
evaluate fidelity easily through prevailing numerical tech- 
niques such as DMRG, ED, and Monte-Carlo simula- 
tions, but also builds a straightforward connection be- 
tween the concepts in quantum information theory and 
those in quantum many-body physics. 
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